BLOCK-ITERATIVE EQUALIZERS 
FOR DIGITAL COMMUNICATION SYSTEMS 



GOVERNMENT SUPPORT 

The Government has rights in this invention pursuant to Contract Number 
N00014-96-1-0930, awarded by the Office of Naval Research. 

FIELD OF THE INVENTION 

The present invention relates to a novel class of equalizers for use in the receivers 
of digital communication systems. 

BACKGROUND OF THE INVENTION 

In the receiver of a contemporary digital communication system, an equalizer 
operates to compensate for various forms of distortion introduced by the propagation 
medium. In wireless communication systems for example, an important source of such 
distortion is time-varying multipath propagation, whereby the transmitted signal travels 
through multiple paths en route to the receiver due to reflections off objects in the 
propagation environment. In wireline communication, such as twisted copper wire pair or 
co-axial cable, such distortion arises due to the frequency response characteristics of the 
physical medium, referred to herein as the "channel", or of other system hardware in the 
network. 

Absent some form of correction or compensation, such distortion can cause 
streams of transmitted bits or symbols to interfere with one another, a phenomenon 
generally referred to as intersymbol interference. In certain code-division multiple-access 
(CDMA) communication systems, distortion may further give rise to a phenomenon 
referred to as inter-chip interference. Intersymbol interference (ISI) is used herein to 
refer to both forms of interference. Precursor ISI is caused by the strictly anticausal 
portion of the equivalent discrete-time channel impulse response, whereby a current 
symbol affects symbols in the past. Postcursor ISI is caused by the stricly causal portion 
of the equivalent discrete-time channel impulse response, whereby a current symbol 
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affects symbols in the future. ISI can lead to severe degradation in system performance, 
especially in wireless settings. Equalizers, designed to compensate for, and to mitigate 
the effects of, such interference, have been under development for more than three 
decades and have evolved to become essential components in virtually all modern 
communication systems. 

It is well known in theory that receivers that realize maximum-likelihood 
sequence detection (MLSD) are asymptotically optimal in bit-error rate at high 
signal-to-noise ratios (SNR), as observed in G. D. Forney, Jr., "Maximum-likelihood 
sequence estimation of digital sequences in the presence of intersymbol interference," 
IEEE Trans. Inform. Theory, vol. IT-18, pp. 363-378, May 1972. However, even when 
implemented efficiently using the Viterbi algorithm, the computational requirements of 
such receivers are prohibitive for even small signal constellations and channels whose ISI 
□ spans only a modest number of symbols. Indeed, although the complexity of MLSD is 

independent of block-length, it does grow as M L ~ l where M is the size of the symbol 
15 alphabet and L is the length of the impulse response of the equivalent discrete-time 
0; channel. As a result, practical equalizers rarely employ MLSD. 

Two approaches have become popular in practice: linear equalizers (LEs) and 
decision-feedback equalizers (DFEs). LEs were developed in the 1960's, and employed 
linear filters to compensate for distortion. Later, in the late 1960's and early 1970's, the 
20 inherently non-linear DFEs were introduced. In both cases, the computational complexity 
is dramatically lower than the theoretically superior MLSD approach - it is essentially 
independent of the size of the symbol alphabet M and proportional to the length of the 
channel impulse response/,. 

Like the LE, the DFE processes the received signal using a linear filter and 
25 makes symbol decisions using a slicer, as disclosed in C. A. Belfiore and J. H. Park, 
Jr., "Decision-feedback equalization," Proc. IEEE, vol. 67, pp. 1143-1156, Aug. 
1979. However, each time a decision is to be made, a linear, strictly causal filter forms 
a weighted linear combination of previous symbol decisions, assumed to be correct, to 
cancel the postcursor ISI at the slicer input. The slicer then generates a decision for the 
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current symbol. Note that once a symbol decision is made, it does not change and is 
used to cancel postcursor ISI when making decisions for future symbols. 

Another equalizer approach, referred to as the "ISI canceler", processes the 
received signal using a linear filter and makes symbol decisions using a slicer, as 
disclosed in Gersho and T. L. Lim, "Adaptive cancellation of intersymbol interference 
for data transmission," Bell Syst. Tech. J., vol. 60, pp. 1997-2021, Nov. 1981. 
However, each time a decision is to be made, a linear noncausal filter forms a weighted 
linear combination of both previous and future tentative symbol decisions made by 
some other equalizer, typically a linear equalizer, to cancel ISI at the slicer input. The 
slicer then generates a final decision for the current symbol. 

LEs are widely used, but often suffer from excessive noise enhancement, even 
when minimum mean-square error (MMSE) design criteria are used. DFEs are capable 
of better performance, particularly at higher SNR, because postcursor ISI is suppressed 
nonlinearly. However, noise enhancement is still an issue in DFEs because precursor ISI 
is only suppressed linearly, and the sequential structure of DFE algorithms make them 
generally incompatible with error-control channel coding due to decoding delay issues, 
despite numerous attempts to merge the two; see, for example, (a) J. G. Proakis, Digital 
Communications. New York, NY: McGraw-Hill, 2nded., 1989. (b) V. M. Eyuboglu, 
"Detection of coded modulation signals on linear, severely distorted channels using 
decision-feedback noise prediction with interleaving," IEEE Trans. Commun., vol. 
COM-36, Apr. 1988. (c) K. Zhou and J. G. Proakis, "Coded reduced-bandwidth QAM 
with decision-feedback equalization," in Proc. Int. Conf. Commun., pp. 12.6.1- 12.6.5, 
June 1988. (d) K. Zhou, J. G. Proakis, and F. Ling, "Decision-feedback equalization of 
fading dispersive channels with trellis-coded modulation," in Proc. Int. Conf. Commun. 
Techn., Nov. 1987. Since error-control coding is required to approach channel capacity, 
an equalizer compatible with this feature is highly desirable. 



SUMMARY OF THE INVENTION 

The present invention represents a major improvement over existing equalizers, 
and its advantages are the result of a fundamentally new equalizer architecture. 

The novel class of equalizers presented herein are configured in a block-iterative 
structure, i.e., they are implemented by a multiple-pass algorithm. An illustrative 
example of the operation of the new class of equalizers is as follows. Let x[n] denote an 
uncoded transmitted symbol stream and r[n] denote a received data stream corresponding 
to the equivalent discrete-time baseband model. On the first pass, this received data is 
processed by a linear filter, and the resulting output is passed on to a slicer, which makes 

a first set of tentative symbol decisions x l [n]. This initial pass is then followed by a 
second pass in which the same received data is processed by a linear filter, and the 
tentative decisions made in the previous pass are used to construct and subtract out an 
estimate of the ISI. The ISI-reduced data is then passed on to the slicer, which makes a 

new second set of tentative decisions x 2 [n] . The process repeats in this manner, until 
adequately reliable decisions are obtained. 

Unlike the conventional DFEs, the novel equalizers of the present invention 
suppress virtually all of the ISI, including precursor and postcursor ISI, in a nonlinear 
manner, significantly enhancing system performance. In fact, theory and simulations 
have demonstrated that the new class of equalizers can achieve a given level of reliability 
(as measured by bit-error rate) using significantly less received signal power than 
conventional equalizers. 

Additionally, unlike the DFE, the equalizers of the present invention can be 
deployed in conjunction with systems employing error-control coding, making them 
attractive for channels operating near capacity. 

A further advantageous feature of the equalizers of the present invention is 
extremely low computational complexity. The order of complexity is comparable to that 
of the DFE, and the adaptive implementation is much less computationally complex. 

In a first embodiment, the present invention is directed to an apparatus and method 
for iteratively equalizing received data transmitted over a data channel in a data 



communication system. At each iteration, a first filter, for example a "feed-forward" 
filter, filters received data according to first filter parameters to generate first-filtered 
data. A combiner modifies the first-filtered data with second-filtered data to generate 
modified data. A decision device, for example a digital slicer generates modified 
tentative decisions based on the modified data, the modified tentative decisions being 
modified with respect to tentative decisions of a previous iteration. A second filter, for 
example a "feed-back" filter, filters tentative decisions from a previous iteration 
according to second filter parameters to generate said second-filtered data. The first 
and second filter parameters are based on the received data, and are preferably 
modified at each iteration. 

In a preferred embodiment, the received data is sampled at a rate higher than a 
symbol rate associated with the received data. The received data may comprise symbol 
data. 

The first and second parameters may be modified at each iteration according to 
channel parameters that are re-estimated at each iteration based on the received data. 

The received data may optionally be encoded for error correction coding, in 
which case the decision device comprises an error-correction decoder, and wherein the 
equalizer further comprises an encoder for encoding the tentative decisions from a 
previous iteration. 

The feed-forward and feed-back filters may comprise linear, non-linear, time- 
variant, time-invariant, infinite- impulse-response (IIR) and finite-impulse-response 
(FIR) types. 

In alternative embodiments, the received data may comprise a plurality of 
received signals received over a plurality of data channels, in which case the equalizer 
further comprises a like plurality of first filters corresponding to the plurality of 
channels. The received data may further comprise combined data for a plurality of 
users, in which case the equalizer further comprises a like plurality of second filters 
corresponding to the plurality of users. 



BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features and advantages of the invention will 
be apparent from the more particular description of preferred embodiments of the 
invention, as illustrated in the accompanying drawings in which like reference 
5 characters refer to the same parts throughout the different views. The drawings are not 
necessarily to scale, emphasis instead being placed upon illustrating the principles of the 
invention. 

FIG. 1 is a graph comparing uncoded quadrature phase-shift keying (QPSK) 
performance predictions for the iterative equalizer of the present invention and various 
10 classical equalizers, based on asymptotic (long ISI) analysis. 

FIG. 2 is a block diagram of a data signal receiver configuration in accordance 
with the present invention. 

FIG. 3 is a block diagram of an equalizer system in accordance with the present 
invention. 



S FIG. 6 is a block diagram of an equalizer embodying the principles of a 

20 multichannel implementation in accordance with the present invention. 
fi FIG. 7 is a block diagram of an equalizer embodying the principles of a multiuser 

O implementation in accordance with the present invention. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
25 The present invention is directed to a novel equalizer system suitable for use in 

modern communication systems. A block-iterative configuration is employed to 
efficiently process received data, in a manner that is reliable and computationally 
inexpensive. 

During each iteration, a set of tentative decisions are made as to the transmitted 
30 data, with the help of tentative decisions from a previous iteration. Reliability of the 
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FIG. 4 is a block diagram of a transmitter that exploits spatial diversity, in 
accordance with the present invention. 

FIG. 5 is a block diagram of a fading channel model for a transmitter array, in 
accordance with the present invention. 
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tentative decisions improves during each iteration, and therefore the tentative decisions 
from the previous iterations become more useful. The iterations continue until further 
iterations do not significantly change the set of tentative decisions. 

Unlike conventional DFEs, the novel equalizers of the present invention suppress 
virtually all of the ISI, including precursor and postcursor ISI, in a nonlinear manner, 
significantly enhancing system performance. In fact, theory and simulations have 
demonstrated that the new class of equalizers can achieve a given level of reliability (as 
measured by bit-error rate) using significantly less received signal power than 
conventional equalizers. In particular, the new class of equalizers is capable of achieving 
a given bit error rate with 2.5 dB less received signal power as compared to a 
conventional system employing a minimum mean-square error DFE on severe ISI 
channels, even for uncoded systems. 

As an example, FIG. 1 compares uncoded quadrature phase-shift keying (QPSK) 
performance predictions for the iterative equalizer of the present invention 20 with 
various classical equalizers, including the minimum mean-square error decision-feedback 
equalizer (MMSE-DFE) 22, the zero-forcing decision-feedback equalizer (ZF-DFE) 24, 
the minimum mean-square error linear equalizer (MMSE-LE) 26, and the zero-forcing 
linear equalizer (ZF-LE) 28, based on asymptotic (long ISI) analysis. The zero-forcing 
(ZF) technique selects the filters that completely eliminate ISI, usually at the expense of 
high noise enhancement. The minimum-mean-square-error (MMSE) technique, on the 
other hand, takes both into account, and attempts to minimize the combined effect of ISI 
and noise enhancement. It is noteworthy that the performance of the iterative equalizer 
of the present invention 20 approaches the additive white Gaussian channel bound, 
asypmtotically achieved by MLSD. When a channel introduces no intersymbol 
interference but only white Gaussian noise, equalization is not necessary. The receiver 
simply uses a symbol-by-symbol decision slicer and the bit-error rate achieved is known 
as the additive white Gaussian channel bound. When a channel introduces both 
intersymbol interference and white Gaussian noise, equalization is necessary. No 
matter what kind of equalization is used, one can never attain a lower bit-error rate than 
the additive white Gaussian channel bound at the same SNR. Receivers that use 



maximum-likelihood sequence detection (MLSD) are asymptotically optimal in the 
sense of minimizing bit-error rate at high SNR over channels that introduce intersymbol 
interference, and so it is expected that the bit-error rate in a receiver using MLSD 
approaches the additive white Gaussian channel bound at high SNR. In FIG. 1, the bit- 
error rate of the iterative equalizer of the present invention 20 approaches the additive 
white Gaussian channel bound at high SNR, which is asymptotically achieved by 
MLSD. 

FIG. 2 is a block diagram of a receiver system 30 for receiving digital pulse 
amplitude modulation (PAM) carrier-modulated data signals transmitted from a far-end 
transmitter (not shown) over a communication channel. 

In the transmitter, a stream of binary data is typically mapped into a sequence x[ri\ 
of complex numbers selected from a discrete symbol set. Multiple bits may be mapped 
into a single data symbol. The sequence of symbols is then combined with the transmit 
pulse shape Ptf(0 to produce the complex-valued continuous-time signal 

*(0 = £*[*]/v('-*rx 

k 

where T is the modulation or baud interval. The real and imaginary components of s(t) 
amplitude-modulate in-phase and quadrature-phase carrier waves respectively, which 
combine to form the passband PAM signal transmitted to the receiver 30. 

At the receiver 30, the received signal Hf) is applied to an input circuit 
comprising several components in series, the order of which may vary from system to 
system. Typically, the input circuit includes a bandpass filter 32 in order to suppress 
energy in the received signal r{t) outside the transmission frequency band. The filtered 
signal 33 is next processed by a phase splitter 34 that rejects negative frequencies and has 
two real-valued outputs representing a complex-valued signal. Alternatively, the 
functionality of the bandpass filter 32 and the phase splitter 34 can be replaced by an 
analytic passband filter, implemented using two physical filters, each with a real-valued 
output. An automatic gain control (AGC) 36 typically follows the passband filtering, 
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which compensates for slow amplitude variations in the received signal caused by the 
channel. A demodulator 38 multiplies the complex-valued signal 37 by 



compensate for the effects of frequency offset and phase jitter. In one embodiment, the 
complex multiply can be implemented using four real multiplies and two real additions. 
Finally, an analog-to-digital (A/D) converter 42 samples the two real-valued signals 39 
representing the complex- valued signal. In high-performance communication systems, the 
incoming analog received signal may by oversampled so that sensitivity to timing and 
synchronization is reduced. A receiver capable of receiving and processing such 
oversampled signals is referred to in the art as a "fractionally-spaced" receiver. Thus, 
samples can be generated at an integer multiple of the baud rate, i.e., at a rate of Q/T 
complex samples/sec, where Q is a positive integer. The baud rate is usually estimated in 
the receiver 30 by a separate timing recovery circuit 48. In this manner, complex 
samples r\n\ are produced, with the index n advancing at a rate of Q/T complex 
samples/sec. The complex samples r[n] can be modelled as 



where a[n] is the equivalent discrete-time baseband model of the PAM communication 
system and where w[n] is additive noise. The demodulator function 38 may optionally 
follow the sampler 42 to allow the demodulator 38 to be implemented in discrete time 



be decision directed, in which case feedback from the output of the receiver is needed. 
The received samples r[n] are next passed to an equalizer system 44 configured in 
accordance with the present invention. 

The equalizer 44 processes the received data r[n] in a block-iterative fashion. 
Specifically, during each iteration or "pass", a linear filter processes the received data, 
and tentative decisions made during a previous iteration are used to construct and subtract 
out an estimate of the ISI. The resulting ISI-reduced data is then passed on to a slicer, 



e K<o c t+0{t)) 



where 9{t) is determined by a separate carrier recovery circuit 40 to 




rather than continuous time. Additionally, both carrier recovery and timing recovery may 



which generates a new set of tentative decisions. With each successive iteration, 
increasingly refined hard decisions (i.e., decisions made by the slicer 58 of FIG. 3 
below) are generated using this strategy. 

FIG. 3 is a block diagram of an iterated-decision equalizer 44 in accordance with 
a preferred embodiment of the present invention. The parameters of all systems and 
signals associated with the / th pass are denoted using the superscript " / " . On the / th 
pass of the equalizer 44, where /= 1,2,3..., the received data r[«] is first processed by 
first filter b l [n] 50, referred to herein as a "feed-forward filter" (preferably in the form 
of a linear filter), producing the sequence: 

v' _ 

r [n] = 2_ Jk b l [kP]r[n-kP] 

where P and Q are preferably prime integers with P < Q. (If a symbol-spaced 
equalizer is to be implemented, then P = Q = 1. If, however, a fractionally-spaced 
equalizer is to be implemented, then P = 1 and Q = 2 in most cases when there is not 
more than 100% excess bandwidth.) The output of a second filter, referred to herein as a 
"feed-back filter" d'[n] 52 and the input to a decision device, for example a slicer 58, 

are preferably both at the baud rate, so r [n] 51 is downsampled at downsampler 54 by 
a factor of Q to provide downsampled signal r'[n] 55. Next, an appropriately 
constructed estimate z l [n] 53 0 f the ISI is subtracted from downsampled signal r'[n] 
55 at subtracter 56 to produce signal x'[n] 57, where 

x l [n] = r l [n]-z'[n] t 

and where 
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Since z l [n] j s intended as an ISI estimate, the zero-delay tap of the feedback filter 
d l [n] 52 is set at zero. The slicer 58 then generates the new set of hard decisions 

A / r i ~~ I r i 

x [n\59 from* [n\ using, for example, the well-known minimum-distance rule. Note 
that alternative functions are equally applicable to the present invention. For example, 
the subtractor 56 may comprise an adder or other means for combining the filtered 
signals from the feed-forward and feed-back filters 54, 52. 

Further details of the invention depend on whether the channel is known at the 
receiver, partially known at the receiver, or unknown at the receiver. The known case 
determines the filter parameters using a fixed set of channel parameters, for example a 
known transfer function of a specific channel (wire, fiber optic, air, etc.). In the 
unknown case, the channel parameters are unknown, and the invention instead determines 
the filter parameters from the received data, using, for example, predetermined training 
symbols. In the partially-known case, an initial estimate of the channel parameters is 
used, and the iterative nature of the equalizer of the present invention is used to improve 
the channel estimate at the end of each pass. In each of the known, unknown and 
partially-known cases, during each pass, the feed-forward and feed-back filter parameters 
50, 52 are modified. 

Known Channel 

When the channel is known at the receiver and P = Q = 1, one possibility is to 
select the feed-forward filter b l [n] 50 and the feed-back filter d l [n] 52 during the / th 
pass to maximize the signal -to- interference + noise ratio (SINR) at the slicer input, in 
which case the discrete-time Fourier transforms of the feed-forward filter b l [n] 50 and 
the feed-back filter d l [n) 52 are respectively 

A* (a>) 

OC 



B l {co) 

n 0 +^o+(^- 1 ) 2 )MH 2 

D l (co) = p'; 1 (A(co)B 1 (cd)-e[a(o))B 1 (co)]) 
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where A{a>) j s the Fourier transform of the equivalent discrete-time baseband model of 
the PAM communication system, s s is the energy of the transmitted symbols, N 0 i s 
the variance of the additive noise due to the channel, E[] is the expectation operator 
(alternatively, the expectation E[-] can be replaced by a frequency average 

< > 2tt J-k md ^ x is me normaliz ed correlation between the transmitted 
data sequence x[n] 59 and the previous set of tentative decisions made on the previous 

iteration * M [«] 60, i.e., Px* * 4**M'* M ["]]/^ . A possible method for 

computing the set of correlation coefficients Px in the special case of M-PSK data 
symbols is the following iterative technique: 



1. Set P x =0 and let / = 1. 

/ 

2. Compute the SINR Y at the slicer 58 input on the / th decoding pass from 



Px via 



C 



Y = 



where 



e x { S )= r u dt 

Js 



and 



e s <\A(o»\ 2 >(l-(p>;>) : 



1-1 
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/ 

Optionally, / can be computed (and in some cases more accurately) from 
Px via 



e s \<A{o))B l (co)> 


2 


N 0 < 


B 1 (<d)\ 2 >+£ s (l-(p'; l ) 2 )<(A(a>)B l (co)-<A(co)B , ((D) >f > 



3. Compute the normalized correlation coefficient P x between the symbols x[n] 
and the decisions x l [n] 59 generated at the slicer 58 via 



where 




4. Increment / and go to step 2. 



Partially Known Channel 

When the channel is partially known at the receiver, the performance of the 
receiver may become particularly sensitive to the accuracy of the channel estimate used 
by the equalizer. To prevent such sensitivity, the equalizer can be modified to use an 
estimate of the channel impulse parameters with some degree of uncertainty as opposed to 
a fixed estimate of the channel parameters (as in the known-channel case). When P = Q 
= 1, one possibility is to use the techniques for the known-channel case described above, 

except that the fixed channel A(co) i s replaced by the channel estimate MjoS) , and the 
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noise variance N 0 i s replaced by N 0 + £ 5 < A(co)- A(o) 2 \ > . Furthermore, it is 

possible to exploit the iterative nature of the equalizer to improve the channel estimate at 
the end of each pass, thus leading to a higher SINR in the slicer 58, which in turn lowers 
the probability of symbol error. At the end of the / th pass, the equalizer has knowledge 

of the received sequence r[n] and the decisions x l [n] 0 f the current pass, for n = 0, 1, 
N-l. For the purpose of estimating the channel impulse response a[n] 9 it may be 

assumed that r[n] and x Y n \ are zero prior to n=0 and after n = N- 1. Furthermore, 
for illustrative purposes, it is assumed here that a[n] is nonzero for n = 0, 1, AM. If 
two vectors are defined 



a / =[a / [0]a / [l]...a / [/-l][ 

x l [n] = [jc 7 ^ 7 ^ - l]...x l [n-L + l]f 

(where the superscript denotes the conjugate-transpose operation for vectors and the 
conjugation operation for scalars, and where the superscript T denotes transposition), 
then the difference between the received data symbols and the response of the channel to 
the most recent set of decisions; i.e., 



00 



r[«]-a 7 x l [n] 



n--cc 



is minimized when 
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where 



00 



0 7 = ^x'^^n] 



n=-oo 



and 



W=-oo 



For M-PSK symbols, the variance of the error in the channel estimate 
e l [n] = a[n] - a l [n] made at the end of the / th pass is approximately 



/ / \ 2 _ 2<|^)| 2 >^(l-^) + N ( 
' *' ~ Ne. 



If, at the / th pass, the variance of the channel estimate made at the end of the 

previous pass is smaller than the variance of the channel estimate used for the (M) th 

pass, then the former shouid be used. If, on the other hand, the variance l ) is 

larger than the variance of the channel estimate used for the (/-l) th pass, then the 
estimate used for the (/-l) th pass is used again for the / th pass. 
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Unknown Channel 

When the channel is unknown at the receiver, an adaptive version of the present 
invention is required, in which optimal FIR filter coefficients are selected automatically 
(from the received data) without explicit knowledge of the channel characteristics. 

The multipass equalizer of the present invention is designed to process received 
data in a block-iterative fashion, so it is ideally suited for packet communication in which 
the packet size is chosen small enough that the channel encountered by each packet 
appears linear time-invariant. As is typically the case with other adaptive equalizers, the 
adaptive iterated-decision equalizer makes use of training symbols sent along in the 
packet with the data symbols. Suppose that a block of white M -ary PSK symbols x[n] for 
/z =0,1,... AM is transmitted; some of the symbols (not necessarily at the head of the 
packet) are for training, while the rest are data symbols. 

The feed-forward filter b l [n] 50 and the feed-back filter d l [n] 52 for the / th 
iteration preferably comprise finite-length filters. Specifically, the feed-forward filter 
b l [n] 50 has J\ strictly anticausal taps and J 2 strictly causal taps plus a center tap, 

while the feed-back filter d*[n] 52 has ^1 strictly anticausal taps and K 2 strictly causal 
taps with no center tap. 

Prior to the first pass ( / = 1), the hard decisions x I/ 7 J are initialized. Since the 
locations and values of the training symbols in x[n] are known at the receiver, the hard 
decisions are initialized as x°[n] = x[n] for the n corresponding to those locations. For 

all the other n between 0 and N - 1 inclusive, x [n\ j s se t to a "neutral" value; as an 
example, for white PSK symbols, this value should be zero. 

On the / th pass of the equalizer where / = 1, 2, 3, ... the slicer 58 

input can be expressed as 
x l [n) = c r q l [n) 
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where 



and 



c' = 



q'M = 



r[n] ■ 



b' = 
d' = 
r[n] = 



b l [-JiP} ■■■ b'l-p) b'[o] b'[p\ ••■ 6'(y 2 p] 



d'[-Ki] •-• d'[-l] d'[l] •■■ d'(A' 3 ] 



r[nQ + JiP] ••• r[nQ + P] r[nQ) r[n<? - P] •■• r[nQ-./ 2 P] 



i'-^n+Ki] ■■■ z'-'fn+l] i'-'ln-l) ±'- l [n - K 3 ] 



Using, for example, a minimum-distance rule, the slicer 58 next generates the hard 

decisions x from x for all n between 0 and N - 1 inclusive, except for those n 
corresponding to the locations of training symbols in x[n] . For those n, the hard 

decisions are set to x l [n] = x[ri] . 

During the / th iteration, there are two sets of data available to the receiver: the 

received data r[n\ and tentative decisions from the previous iteration X [n], 

n=0,\,..., N-l. If it is assumed that x[n] » x ~ [n] 

for the purpose of determining 
the optimal filters (as is similarly done in the adaptive DFE in decision-directed mode), 

then it is reasonable to choose the feed-forward filter parameters b l [n] 50 and the feed- 
back filter parameters d l [n] 52 so as to minimize the sum of error squares: 
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00 



n=-co 



Since this is a linear least-squares estimation problem, the optimum C is represented 
by: 



where 



00 



« = -00 



and 



»' = Zi'-t-Ji'W. 



«=-oo 
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The resulting equalizer lends itself readily to practical implementation, even for large 
filter lengths. In particular, the matrix <J>' can be efficiently computed using correlation 
functions involving r[n] and x l [n] . Specifically, if the following deterministic correlation 
functions are defined as: 



P[«l,"2] 
if/'[«] 



£ r[kQ- ni ]r*[kQ- n2 \ 

k - -oo 
k = -oo 
k - -oo 



where 



x l -\nlQl n = 0,Q,2Q,...,(N-\)Q 
0, otherwise 



then O 7 can be expressed as 



P V 1 
W 1 * Q l 



where 



P = 



ftjf, -Jf\ 



PW 



P[ -Jf, Jf] 



19 



O'[0] 
9>[-Q] 



0'[Q] 
O'{0) 



0' = e'HAVUQ] e'(-(A',-2)Q] 

6'[-{Kx+l)Q] 0'[-K,Q\ 
e>{-(Ki+!<iiQ} 6' [-(A'i+A' 5 -1)Q) 



•• 0'[(A'i-2)Q] 

fl'(0] B'[2Q] 
8'[-7Q\ e'[o] 



e'((A',+/< 2 )Q] 
fl'[(A',-=-A' 2 -l)Q] 

9'[(A' 3 +1)Q] 
5'[(A' 2 -1)Q] 

e'[o) 



and 



0'[_./ 2 p_A' 1 Q) v^' (— ^y 3 — (A'i— lyQ] ■■• ^'[-J^P-Q] v'[-7;P-K?] ••• i>\-JiP^-K-C\ 



Note that in general, none of the matrices P, ©' , or i s Toeplitz. However, ©' is a 
(K l +K 2 + 1) x (^j +tf 2 + 1) Toeplitz matrix with its (K { + 1) th row and (K l + 1) th column 
deleted, and P exhibits cyclostationary properties with each diagonal having a period of Q. 
Moreover, depending on the values of P and Q, the matrix *F 7 may have a significant 

number of entries that have the same value. Thus, the special structure in can be 
exploited so that the redundant computation of repeated entries is avoided. Furthermore, since 
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is a Hermitian matrix, it is only necessary to compute half of the distinct entries. In the 
special case of p = Q = i, P is a Toeplitz matrix, and V F / can be thought of as originating 
from a (J\ +^2+ l ) x (^l +^2 + l ) matrix with the elements on each diagonal being equal, but 

now with the (K Y + 1) th column deleted. In this special case, the computation of requires 
only 2(y } + J 2 +K X +K 2 + 1) of the ( J x + J 2 + K { + K 2 + 1) 2 entries, as there are only that many 
distinct entries. If the total number of taps in the equalizer is 100, then the computation of 

has been made approximately 50 times more efficient! The computation of 
also be made more efficient. Specifically, since can be expressed as a partitioned matrix in 
which P and & l are square, can be partitioned into 



can 



K]" 1 = 



pi q,i 



where P' , & , and have the same sizes as P, and respectively. Note that this 
takes advantage of the fact that the inverse of a Hermitian matrix is also Hermitian. Then, the 
matrices p , and T 7 can be computed either by the formulas 



©' = [©'j- l + (*'[0']- 1 )t-(/3-*'[0']- 1 *")-'. (*'[©']-!) 



or by the equivalent formulas 
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P l = ^" 1 + (/3- 1 ^)-(©'-* /, /3- 1 *')- 1 .(/3- 1 ^')f 

(See H. Lutkepohl, Handbook of Matrices. Chichester, England: Wiley, 1996.) The factors 
within parentheses in the two above sets of equations indicate quantities that need to be 
computed only once for a given /. So if the first set of formulas is used, then 



^[0'j- 1 , (P - *'[0 , ]- 1 *' , )- 1 > and 00 - *'[0']-^")-i . (^[0'J- 1 ) 

should be computed in that order during each iteration. Similarly, if the second set is used, 
then 

p-\ (/r 1 *'), (©' - ^'/r 1 *')- 1 , and ■ (©' - #"/r 

should be computed in that order during each iteration. In either case, the task of computing 




has been modified into a task involving the computation of an inverse matrix the size 



of P and another inverse matrix the size of ©' . One advantage of using the second set of 
formulas is that the same /T 1 is used during all iterations, so it only needs to be computed 
once during the whole multipass process. Thus, the only matrix inverse that needs to be 
computed during each iteration (except for the first iteration) is 
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• 



and the original task of computing the inverse of a 

(J\ +J 2 +K { +K 2 + 1) x (Jj +J 2 + K X +K 2 + 1) matrix has been converted into a task involving 
the computation of the inverse of a (AT, x (K x +^2) matrix. If «A is large, and 

J x = J 2 = K x = K 2 , then a reduction of computation by an order of magnitude is plausible. 
We now turn to a couple of implementation issues. First, when P = Q = 1 , the finite-length 
adaptive filters ideally approximate 



n'r \ A* (a) 

B (a)) oc ^ 

n 0 +^(i+(p^)MM 2 

D l {co) = P 1 ; x (a{(d)B 1 (cq) - e[a(co)B 1 (co)]) 



which are infinite in length. The optimal feed-forward filter b'[n] includes a filter matched to 
a[n] , and the optimal d l [n] includes a cascade of a[n] and the corresponding matched filter, 
suggesting that a reasonable rule of thumb is to select J X ,J 2 ,K X ,K 2 such that each FIR filter 
covers twice the span of the channel impulse response. Second, the block- iterative nature of 
the equalizer allows the training symbols to be located anywhere in the packet. Since (in 
contrast to the DFE ) the locations do not appear to affect equalizer performance, the training 
symbols are arbitrarily chosen to have uniform spacing within each packet. 



Error-Control Coding 

In most communication systems, it is preferred to have the ability to apply error-control 
coding (ECC) to the data prior to modulation and transmission. Examples of error-control 
codes include block codes, convolutional codes, and trellis-coded modulation. The present 
invention includes two implementations of the equalizer that allow such coding to be explicitly 
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taken into account at the receiver. In one version, equalization and error-correction decoding 
are performed separately. This implementation is obtained by processing the final decisions of 
the equalizer using a decoder placed after the equalizer in FIG. 2. In the other version, 
equalization and error-correction decoding are performed jointly, leading to much improved 
5 performance. This implementation is obtained by generalizing the iterative equalizer of FIG. 3 
by coding the symbols x l [n] using an error-correction encoder before being filtered in the 
feedback loop by d l [n] , and by replacing the slicer depicted in FIG. 3 with an appropriate 
soft- or hard-decision error-correction decoder. For example, if trellis coding (see (a) G. 
Ungerboeck, "Channel coding with multilevel/phase signals," IEEE Trans. Inform. Theory, 
io vol. IT-28, pp. 55-67, Jan. 1982. (b) G. Ungerboeck, "Trellis-coded modulation with 

redundant signal sets Part I: Introduction," IEEE Commun. Mag., vol. 25, pp. 5-11, Feb. 
1987. (c) G. Ungerboeck, "Trellis-coded modulation with redundant signal sets Part II: State 
0 of the art," IEEE Commun. Mag., vol. 25, pp. 12-21, Feb. 1987.) were employed, a 
[jj classical soft-decision Viterbi decoder could be used in the resulting receiver. As another 
[If 15 example, if turbo-type codes (see (a) C. Berrou, A. Glavieux, and P. Thitimajshima, "Near 
|j Shannon limit error-correcting coding and decoding: Turbo codes," in Proc. Int. Conf. 
Jjj Commun., pp. 1064-1070, 1993. (b) C. Berrou and A. Glavieux, "Near optimum error 

, correcting coding and decoding: turbo codes," IEEE Trans. Commun., vol. 44, pp. 

O 

y 1261-1271, Oct. 1996.) were used in conjunction with this equalizer, a doubly- iterative 

Fjfio composite receiver could be used in which iterative turbo-decoding takes place in lieu of the 

□ slicer as a nested subsystem of the larger iterative equalizer structure. 

Implementation for Transmitter Antenna Arrays 

Many communication systems, especially wireless ones, use a transmitter antenna 
25 array, whereby the data to be transmitted are sent from each of a set of different antenna 

elements. One possible transmitter structure which uses m antennas to exploit spatial diversity 
is depicted in FIG. 4 (see G. W. Wornell and M. D. Trott, "Efficient signal processing 
techniques for exploiting transmit antenna diversity on fading channels," IEEE Trans. Signal 
Process., vol. 45, pp. 191-205, Jan. 1997). For convenience the focus of the present 
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discussion is on equivalent discrete-time baseband models. The symbol stream x[n) 72 is 
scaled down by a factor of \Tm at scaler 74 and then processed by a linear precoder 76 at each 
of the constituent antennas of the transmitter to produce the transmitted symbol streams 

v \[n], v M, v,„[«] 78. 

The present invention includes an implementation of the block-iterative equalizer that is 
designed for use with a transmitter antenna array like the one described above. If the channel 
can be described by a nonselective fading model, shown in FIG. 5 , then the signal at the 
receiver 82 can be modelled as 

m 

r[n] = £ OjVjln] +w[n] 
i = l 

where w[n] 84 is additive noise and a h a 2 , a,„ 80 are fading coefficients. Alternatively, 
the received signal can be expressed in the form 

r[n] =^x[k]a[n-k] +w[n] 
k 

where 

1 m 

a[n] = —Ya i h i [n] 
ylm i = l 

is the impulse response of the "effective" channel generated by the antenna precoder. The 
channel has a frequency response of 

1 m 

A(o) =— £ 

so the antenna precoding effectively transforms the original nonselective fading channel into a 
frequency-selective fading channel, and the implementation of the block-iterative equalizer 
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corresponding to FIG. 3 and described above is applicable. As discussed above, specific details 
of the implementation depend on whether the "effective" channel is known, partially known, 
or unknown. 

Multichannel Implementation 

In many communication systems, especially wireless ones, multiple copies of the 
transmitted waveform are received, each being distorted differently. Such is the case, for 
example, when the transmission is received by a multiple element array. In such cases, it is 
common to use an equalizer with what is referred to as a multichannel implementation for 
jointly processing the received streams. For such applications, the present invention includes 
an efficient multichannel implementation of the block-iterative equalizer, which can enhance 
performance over the single channel implementation. 

The structure of a multichannel, fractionally spaced iterated-decision equalizer is 
depicted in FIG. 6. For each of the m channels, the equalizer receives a signal that is typically 
modelled as 

r M = ~k] + Wj[n] for i = 1, 2, m 

k 

where x[n] is the sequence of transmitted symbols, a,[«] is the equivalent discrete-time 

baseband model of the i th channel, and is the additive noise of the / th channel. 

Therefore, in this embodiment, the main difference in the equalizer is that instead of one 
fractionally-spaced feed-forward filter, there are now m such filters 90. The operation of the 
equalizer remains essentially the same, except that the output of the feedback branch 92 is now 
subtracted from the sum of the feedforward branch outputs, rather than from just a single 
feedforward branch output. 

If the multiple channels are unknown at the receiver, the technique employed in the 

single channel case can be generalized. Thus, the slicer input x l [n] can be expressed in 
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vector form as 



5 where 



x l [n] = c l <\ l [n\ 



10 



O 



01 



0J15 



03 



Q 

w 
m 



** 2 0 



and 



c = 



-d ! 



q'[n] = 



r iN 

r 2 [n] 

r m [n] 
x'- l [n] 



bi = 




6((-P] 6{[0] b<[P\ ••• 


6;.(J 2 P1 










d' = 


d'i-Kt) ••• 


<f'(-l] ■•• d'{K,\ 




r,-(n] = 


r,(n<?-r./ l P] 


■ ■■ r;{nQ+P] Ti[nQ] r, 


[«Q - P) 


x'" l (n] = 


i'- l [n+ A',) 


... i'-»(n+l] £'-'["- 1) •• *' 



1* 



If A W.A5W,. ..A>], and</'[#i] are chosen so as to minimize the sum of error squares: 



00 



s(c')= £ S'-'M-c'Ym 



n=-oo 



25 then the optimum is 
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where 

00 

and 

00 

u' = E* M W["]. 



The computation of the matrix O' can be made more efficient by first defining the following 
deterministic correlation functions 

oo 

Pij[n u n 2 ] = E n[kQ - n x )T][kQ - n 2 ] 

k—-oo 

CO 

*'N = E y l - l m- v [k-n) 

h=~co 
oo 

k=-co 

where 



0, otherwise 




y l ~ l [n] = \ 
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Thus , for m channels, the matrix o' may be rewritten as 



011 012 

(3\ 2 P22 



01m *i 

/3 2m *7 



where 



fl'[0] 

o'l-Ql 



o'[o] 



fl'[(/fi-i)Q] o'[(/m-m)Q] 



0' = 



fl'l-(ffi+i)Q] 



0 ! [0] 

0'(- 2 Q] 



0' [20] 



9' 10] 



and 



4(JiP-Q] ^!t-/iP^?l 
^'.[(J,-l)P-Q] ^SlC/i-iJP-^l 



^SlJiP+KsQl 

t/aCA-iy+K*?) 
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Note that in general, none of the matrices Py, or } s Toeplitz. However, ©'is a 
(£, +K 2 + 1) x (K l +K 2 + 1) Toeplitz matrix with its + 1) th row and (K { + 1) th column 
deleted, and Pij exhibits cyclostationary properties with each diagonal having a period of Q. 
Moreover, depending on the values of P and Q, the matrix may have a significant 

number of entries that have the same value. Thus, the special structure in O 7 can be 
exploited so that the redundant computation of repeated entries is avoided. Furthermore, since 

is a Hermitian matrix, it is only necessary to compute half of the distinct entries. 

As in the single-channel case, the computation of [<&'] for the multichannel equalizer 

can be made more efficient by partitioning O 7 into four submatrices such that the upper left 
and lower right submatrices are square, and then applying the formulas for the inverse of a 
partitioned matrix. Since these formulas involve computing the inverses of two matrices with 
the same sizes as the upper left and lower right submatrices, it is usually best to partition 

$ into submatrices of approximately the same size so that computation of a large matrix 
inverse is avoided. 



Multiuser Implementation 

The present invention further comprises a block- iterative multiuser equalizer that jointly 
takes into account multiple-access, co-channel, and extra-network interference, which are 
generally experienced in communication networks, especially wireless networks. For 
multiple-access interference in particular, the invention jointly suppresses intersymbol and 
multiple-access interference in a manner generalizing the structure of FIG. 3. 

The signal at the receiver is typically modelled as 

P 

/ = ! k 
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where p is the number of users, x t [n] is the sequence of symbols transmitted by the # th 
user, a^n] is the equivalent discrete-time baseband model of the channel encountered by the 
'th user, and w[n] represents additive noise and interference. 

The multiuser, fractionally spaced iterated-decision equalizer consists of p main 

branches, with each branch designed to decode the symbol stream of a different user. The 
structure of the / th branch, shown in FIG. 7, is similar to the single-user equalizer shown in 
FIG. 3 except that the present embodiment includes an additional p - 1 feedback loops 94 to 
cancel multiple-access interference caused by the other p - 1 users. Because each branch 
requires access to the tentative decisions made in the previous iteration for all p users (i.e., 

each branch requires access to x{~ \n] ,x{~ \n] ,...,x p ~ \n] ), all branches must complete 
the (/- 1) th iteration before any branch proceeds to the / th iteration. 

If the channels are unknown at the receiver, the approach for the single-user case can 

be generalized for the multiple-user case. Thus, for the / th user, the slicer input */[«] can be 
expressed in vector form as 

*/["] = cfq,'[«] 

where 





b{ 




r[n 


] 








x'f 1 






-4 


q!W = 




[n] 








_*ir l 
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and 



10 



0 
4i 

SI 



< 



r(n] = 



b[[-JiP) ■■ 
r[nQ -f ii P] 



•• r(nQ + P] r[nQ] r[nQ - P] ••■ r(nQ-J 2 P] 



if j = i 



_ ij-'ln + ^i] •■ i'- l [n-l] ... x;.-'(n-K 2 ] J , 



if ;' = x 
if;#i 



where; = 1,2,...,/;. If Z>/[/7] ,d l u [n] ,d l 2j [n] , . ..J pi \ ri\ are chosen so as to minimize the 



£320 sum of error squares: 



00 



*/~ ! m-c; q;w 



n=-oo 
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then the optimum c / for the i th user i 



is 



cU=[of]"'u; 
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where 

00 

^!=Xq![«]qf[«] 

n=-ao 

and 

M=-00 

Joint Transmitter and Receiver Design 

The present invention is further directed to an embodiment in which the transmitter has 
available at least partial information about the channel state (for example, via a low-delay 
feedback path from receiver to transmitter). In this situation, significant performance 
enhancements can be achieved. For example, if the transmitter has exact knowledge of the 
channel characteristics, part of the equalization can be performed at the transmitter so that 
noise enhancement and/or error propagation can be avoided. This idea is analogous to 
Tomlinson-Harashima precoding; see (a) M. Tomlinson, "New automatic equalizer employing 
modulo arithmetic," Electon. Lett., vol. 7, pp. 138-139, Mar. 1971. (b) H. Miyakawaand 
H. Harashima, "A method of code conversion for a digital communication channel with 
intersymbol interference, " Trans. Inst. Electron. Commun. Eng. Japan, vol 52-A, 
pp. 272-273, June 1969. (c) H. Harashima and H. Miyakawa, "Matched-transmission 
technique for channel with intersymbol interference," IEEE Trans. Commun., vol. COM-20, 
pp. 774-780, Aug. 1972. 

While this invention has been particularly shown and described with references to 
preferred embodiments thereof, it will be understood by those skilled in the art that various 
changes in form and details may be made therein without departing from the spirit and scope 
of the invention as defined by the appended claims. 

For example the first and second filters specified above may comprise linear, non- 
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linear, time- variant, time-invariant, infinite-impulse response, and finite impulse response 
filters. Although the above description of the preferred embodiment of the present invention 
pertains to equalization of digital signals, the principles of the present invention apply equally 
well to equalization of analog signals, using, for example, analog equivalents of the first and 
second filters, the combiner, and the slicer. 
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